Atomic homodyne detection of continuous variable entangled twin-atom states. 
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Historically, the completeness of quantum theory has been questioned using the concept of bipar- 
tite continuous variable entanglement [T]. The non-classical correlations (entanglement) between the 
two subsystems imply that the observables of one subsystem are determined by the measurement 
choice on the other, regardless of their distance. Nowadays, continuous variable entanglement is 
regarded as an essential resource allowing for quantum enhanced measurement resolution [2], the 
realization of quantum teleportation 3-5 and quantum memories [3j [6] , or the demonstration of 
the Einstein- Podolsky- Rosen paradox [UEHS]. These applications rely on techniques to manipulate 
and detect coherences of quantum fields, the quadratures. While in optics coherent homodyne de- 
tection [10] of quadratures is a standard technique, for massive particles a corresponding method 
was missing. Here we report on the realization of an atomic analog to homodyne detection for the 
measurement of matter- wave quadratures. The application of this technique to a quantum state 
produced by spin-changing collisions in a Bose-Einstein condensate [TT] [12] reveals continuous vari- 
able entanglement, as well as the twin-atom character of the state [13] . With that we present a 
new system in which continuous variable entanglement of massive particles is demonstrated [6] [TJ . 
The direct detection of atomic quadratures has applications not only in experimental quantum atom 
optics but also for the measurement of fields in many-body systems of massive particles [15] . 



Continuous variable entangled states, whose inter- 
mode entanglement is reflected in quadrature as well as 
population correlations, are routinely generated by para- 
metric downconversion in optical parametric amplifiers 
(OPA) [10]. Spin-changing collisions in ultracold bosonic 
quantum gases provide a similar nonlinear process for 
matter-waves generating atom pairs with spin up and 
down - twin-atoms [T6HT9] . Correlated photon or atom 
pairs in the signal | t) and idler | \) modes are created 
from a pump field via nonlinear interactions described 
by the Hamiltonian hao(a^a\. + a^a^), where a\ is the 
creation operator of the respective mode k and 2tyH is 
Planck's constant. This Hamiltonian presumes a large 
amplitude coherent pump field, whose depletion or distor- 
tion by the signal and idler modes is negligible such that 
it can be treated classically. In optics the effective non- 
linearity ao is proportional to the amplitude of this pump 
field and the susceptibility of the medium. In the analo- 
gous regime for atoms the effective nonlinearity originates 
from interatomic interactions and is proportional to the 
population in the pump mode, that is, the Bose-Einstein 
condensate (BEC). For initially empty signal and idler 
modes quantum fluctuations are amplified and the out- 
put is the so called two- mode squeezed vacuum [TO] , This 
state is characterized by the vanishing of the mean field 
amplitude (Ek) = in each of the modes, while the mean 
intensity is nonzero (Ik) oc (E^Ek) > 0. The associated 
two-mode entanglement is revealed by the quadratures 
X±((p) = a±e~ lip ± ale 1 ^ + h.c, the signature being the 
squeezing of their variance for appropriately chosen ip. 
To access the two-mode quadratures the coherent ho- 
modyne measurement technique is employed, whereby 



the conjugate canonical quadratures Xk = a\ J r ak and 
Yk = i(a\ — dk) are measured by mixing the quantum 
mode with a strong classical reference field, the local os- 
cillator a^o ~ VNhO e ~ l<p with large amplitude V^Vlo- 
Experimental control of the local oscillator phase ip pro- 
vides access to the continuous distribution of single mode 
quadratures Xk(cp) = e %ip a\ + e~ l(fi dk and their two- mode 
counterparts X±(tp) (Xk and Yk are used as abbrevia- 
tions for Xk(0) and Xk(^)). For the squeezed vacuum 
state the quadrature distribution of individual modes is 
isotropic and centered around the origin reflecting the 
undefined phase. The correlation of the modes leads to 
anisotropic and squeezed distributions of the two-mode 

quadratures X±(ip) as seen in figure [tjl ma edi . 

Here we report on the development of atomic ho- 
modyne detection and apply this technique to measure 
quadratures of twin-atom matter-wave fields generated 
by controlled spin-changing collisions in a BEC. The ex- 
periments involve a spinor BEC of Rubidium 87 trapped 
in a ID optical lattice with few hundred atoms per site 
and a lattice spacing of 5.5 /im (Fig. [TJd). A high lattice 
potential assures that tunneling between sites is negli- 
gible on the experimental timescale such that the con- 
densates in the different sites are independent [21 . The 
density in the center of the lattice sites is in the or- 
der of 2 x 10 14 cm -3 resulting in a minimal spin healing 
length of £ ~ 1 /im, which is comparable to the exten- 
sion of the on-site wave function (approximately 1.1 /im 
FWHM). The spatial degrees of freedom are therefore 
frozen and the dynamics happens exclusively in the hy- 
perfine spin [22]. The condensate is initially prepared 
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FIG. 1. Analogy to optics and measured population correlations of twin-atom states, a, Parametric downconver- 
sion with light and the analogy to atomic spin changing collisions. In a nonlinear medium effective interactions result in the 
creation of photon pairs in the signal (red) and idler (blue) modes from pump mode photons (green). The quadratures of the 
output modes are centered around the origin - the individual ones are isotropic, while the two-mode quadratures are squeezed 
reflecting their correlations (purple). Pair creation due to spin changing collisions in tight traps is an analogous process in 
quantum atom optics, b, Experimental system and exemplary raw absorption pictures. The illustration on the left shows 
schematically the ID optical lattice which is superimposed by an optical dipole trap to tightly confine the atoms in all direc- 
tions. In each site selected Zeeman levels of the hyperfine spin are the only accessible degrees of freedom. The states involved 
in the spin changing collision process are highlighted. Three typical experimental pictures after spin changing collisions are 
presented on the right where the clear spatial separation of sites and states can be seen, c, Population correlations. Population 
correlations are visible directly on the raw data where the atom numbers in the signal and idler modes fluctuate strongly (blue, 
red) while their difference (black) shows small noise. The left part of the graph is a zoom into 25 experimental realizations 
clearly showing that the fluctuations are common mode. 



in the (F, ra^) = (2,0) hyperfine state serving as the 
analog to the pump mode in an OPA. Making use of 
a combination of quadratic Zeeman shift and state de- 
pendent microwave dressing [23] [24] the spin dynamics 
is energetically restricted to a spin 1 subspace defined 
by the 0, ±1 Zeeman states (see supplementary informa- 
tion). We choose a constant evolution time under spin 
changing collisions of 22 ms in all our experiments. In 
the regime of small pump mode depletion spin-changing 
collisions can be described by an OPA type pair coupling 
between the Zeeman states [TTJ [12] [T6l [17] with the +1 
and —1 states as signal (t) and idler (|) modes. In op- 
tics the integrated nonlinearity is restricted by short in- 
teraction times, a limit which is surpassed in the atomic 
system, such that continuous variable entangled states 
with large population can be realized. Since the under- 
lying process generates pairs of atoms in the signal and 
idler modes - twin-atoms - the mean and variance of 
the population difference 7V_ = N± — should ideally 
vanish although the total population N+ = N± + fluc- 
tuates strongly. This can be experimentally confirmed 
by absorption imaging of the atomic cloud [21] where 
state selectivity is achieved by Stern-Gerlach separation 
(Fig#, c). 

A quantitative analysis of the variance of the popula- 
tion difference A 2 7V_ and the distribution of the popula- 



tion sum N+ gives a first indication of the quantum state 
of the system (Fig. [2]). Due to the pump mode population 
dependent nonlinearity the mean atom number in signal 
and idler mode (iV+) increases nonlinear ly with the total 
atom number iV, resulting in a growing fraction (iV+) /TV. 
In the small N < 300 limit the observed distribution of 
the total population N+ matches the prediction for the 
squeezed vacuum (taking detection noise into account by 
convolution with a gaussian) with maximal squeezing pa- 
rameter r ~ 2 [10]. The corresponding mean population 
is (N+) « 25 and the data shows a comparably large stan- 
dard deviation - in contrast, the noise in the population 
difference 7V_ almost vanishes (AN- < 1.9 atoms). It 
is important to note that the atomic variances reported 
throughout this manuscript are corrected for detection 
noise as detailed in the supplementary information. For 
larger total atom numbers TV the pump depletion can- 
not be neglected any more resulting in a breakdown of 
the analogy to the optical parametric amplification pic- 
ture and the observed distributions of iV+ differ from the 
squeezed vacuum (insets in Fig. [2]) [22]. In this regime 
the well controllable spin dynamics offer prospects for 
the deterministic generation of non-gaussian entangled 
twin-atom states [25 . Density dependent two-body spin 
relaxation loss deteriorates the perfect pair correlations 
such that the relative population noise grows with N. 
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FIG. 2. Population of the signal and idler modes. The red data shows the mean population (iV+) = (N^ + N±) 

versus total atom number N. The bars indicate the fluctuations (1-s.d.) of the data. The large uncertainty is a feature of the 
nonlinear process itself which leads to a non-gaussian distribution of the mode population. The gray area is the sub-poisson 
regime and its upper bound is obtained by a spline fit to (N+). Blue data points show the variance of the population difference 
A 2 iV_, which is compatible with zero for small N (bars indicate the 1-s.d. statistical uncertainty). Its increase with N is 
reproduced when taking particle loss due to spin relaxation into account (black line). The insets show the distribution of the 
mode population for indicated total atom numbers. For small N the distributions match the prediction for a non-depleted 
pump with the measured mean population (iV+) (black lines) while for larger N they clearly differ. The fitted squeezed vacuum 
distribution for 250 < N < 300 corresponds to a squeezing parameter of r « 2 [10] . 



Averaging the results for all total atom numbers we find 
the noise A 2 7V_ suppressed by 6.9 dB below the poisson 
level (7V+). 

The observable populations do not suffice to charac- 
terize the two-mode state; we must access coherences be- 
tween the modes, the two-mode quadratures. We em- 
ploy atomic homodyning by analogy to the experimental 
technique used in optics, whereby the quadratures are 
measured relative to a reference field. Yet in the atomic 
case the population of available local oscillator (LO) ref- 
erence fields is often limited to rather small atom num- 
bers. This has important consequences for the observ- 
able quadratures X k = (a^a L o + a k^ LO )/yJ («LO a Lo), 
which in general differ from the canonical ones. A mea- 
sure for this difference is given by the ratio of the quan- 
tum mode population to the local oscillator population, 
byjvhich the Heisenberg uncertainty relation is modified 
AX k AY k > 1 - (N k )/(al a L0 ) [26]. 
For the twin-atom state one expects peculiar coherence 
properties: The quadratures of the individual modes fluc- 
tuate strongly around a vanishing mean, while the two- 
mode quadratures X±(ip) show reduced fluctuations for 
certain ip reflecting the quantum correlations. To re- 



veal these characteristics in the coherences of the gen- 
erated quantum fields we employ a homodyning scheme 
such that information about the two-mode quadratures 
X±(ip) can be obtained. Measurement of these observ- 
ables requires comparison of each mode with a local os- 
cillator. We ensure relative local oscillator phase sta- 
bility by using the pump mode (2,0) as the simultane- 
ous reference, whose phase ip can be controlled by mi- 
crowave dressing. An atomic three-port beam splitter is 
realized by radio-frequency coupling of signal and idler 
modes to the pump. Analysis of the fluctuations in atom 
number sum and difference N^,-i) ^ ^(2,1) between the 
(2, —1) and (2, 1) Zeeman states after the mixing pro- 
vides, after proper normalization, an upper bound (ub) 
for the variance of the two-mode quadratures A 2 X± b ((/?). 
The measured fluctuations include further noise contri- 
butions governed by population fluctuations in the signal 
and idler modes, which are small only for the quadrature 
difference. Further details about the analysis in the three 
mode framework and the choice of the normalizations can 
be found in the supplementary information. Figures [3^i 
and b show raw data in the regime of small pump deple- 
tion of the measured normalized atom number difference 
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FIG. 3. Atomic homodyning. a, Homodyne measurement of the quadrature difference. After symmetric coupling to the 
wif — state (inset in b) the normalized population difference in the = ±1 states provides a measure for the two-mode 
quadrature difference. Phase dependent noise is visible directly on the raw data while no dependence is observed in the mean 
population difference (solid blue line), b, Homodyne measurement of the quadrature sum. An upper bound for the variance 
of the two-mode quadrature sum is obtained from the fluctuations of the normalized population sum after the coupling. The 
fluctuations of the raw data are strongly phase dependent while the mean (blue solid line) stays constant. The variance minimum 
is shifted by ^ as compared to a indicating noise suppression in the orthogonal two-mode quadrature. The data shown in a and 
b corresponds to all measurements with total numbers 150 < N < 200 and the red lines indicate the inferred standard deviation 
above and below the mean after subtraction of the readout noise, c, Classical vs. quantum preparation - inter-site correlation 
measurement. Characterization of classical phase noise is crucial to obtain information about the individual quadratures of the 
signal and idler modes. For initial states prepared classically that is by linear coupling, strong inter-site correlations between 
different sites in the optical lattice after mixing with the local oscillator (illustrations) reveal classically washed out coherences 
(top panel). In contrast, when the (F = 2,?tif = ±1) states are populated by spin changing collisions no correlations are 
observed within the statistical uncertainty showing that (X^)) = as a direct result of the spin dynamics. 



and sum versus local oscillator phase ip. No phase depen- 
dence is visible in the mean, nevertheless, the fluctuations 
are strongly modulated indicating phase correlations be- 
tween the two atomic modes. 

We set up a different experiment to obtain the quadra- 
ture fluctuations of the signal and idler modes individu- 
ally. As local oscillators we prepare population in the 
(1, =Fl) states before initiating the dynamics. Mixing 
of the local oscillators with the quantum modes is done 
by a two-photon microwave and radio-frequency | cou- 
pling pulse. For both, the signal and idler mode, we find 
coupling phase independent fluctuations of the atomic 
quadratures X^-m = (iV(2,±i) — ^(i,=fi))/\/ (^lo) an d a 
vanishing mean. The mean local oscillator population 
(^V"lo) °f the (1, =pl) state was measured by omitting the 
mixing pulse. 

To assure that our observations are indeed resulting di- 
rectly from the spin changing collision process the tech- 
nical phase noise has to be characterized. The ID op- 
tical lattice setup allows for this characterization by a 
reference measurement in which the population in the 
(2, ±1) states is also prepared by electromagnetic cou- 
pling and spin changing collisions are tuned off-resonance 
during the evolution time. The inter-site correlations ob- 
served after the mixing are shown in figure [3J3 for this 
reference experiment (top) and for the quantum state 



prepared by spin changing collisions (bottom). In both 
cases we observe similarly large on-site fluctuations in the 
atom number difference. However, in the reference ex- 
periment we can use the inter-site correlations to reduce 
these fluctuations by a factor of 20 such that the remain- 
ing fluctuations are close to the expected noise limit for 
two coherent modes. This shows that inter-site correla- 
tions detect finite coherences that have been washed out 
by classical phase noise in our measurements. In case 
of spin changing collisions no correlations are observed 
within the statistical uncertainty, meaning that the van- 
ishing coherences are indeed caused by the process it- 
self. Given this result, technical phase noise does not 
influence the sum of the orthogonal quadrature variances 
A 2 X^(|) + A 2 Yf(|) = 2A 2 X^(|) (the measured quadra- 
ture variance is LO phase independent), such that it is 
a useful observable to characterize the quadrature fluc- 
tuations. For further details on the inter-site correlation 
measurement we refer to the supplementary information. 

Inter-mode entanglement of the twin-atom state is re- 
flected in suppressed noise in the two-mode quadratures 
compared to the quadrat urejluctuations of t4ie individual 
modes A 2 X_+A 2 Y+ < A 2 X t +A 2 X ; +A 2 F t +A 2 y i [13]. 
When investigating this inequality in our experiment, 
care has to be taken to assure comparability of the single- 
and two- mode measurements [26] . In both measurements 
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FIG. 4. Two-mode quadrature fluctuations and mode inseparability, a, Two-mode sum and difference quadratures. 
The measured variances A 2 X+ b ((^) (blue) and A 2 X^°(ip) (red) - upper bounds for the sum and difference two-mode quadrature 
variances - are plotted versus local oscillator phase for total atom numbers 150 < N < 200. The shown experimental results 
are not corrected for readout noise (in contrast to the values stated throughout the text and in b, c) but its magnitude is shown 
as the black line, b, Continuous variable quadrature entanglement. Inseparability is detected in the gray shaded region where 
EA 2 ., the upper bound for the sum of the two-mode quadrature variance minima (green), is smaller than EA 2 , the sum of the 
variance in the individual quadratures (black). The inset shows the ratio EA|/EA? . c, Minimum of the two-mode difference 
quadrature fluctuations. For small total atom numbers N - in the non-depleted pump regime and where spin relaxation loss 
is small - the minimum two-mode quadrature variance A 2 (ipo) inferred from a local quadratic fit is at the noise limit for 
two coherent fields (dashed black line). The inset shows a zoom into the small N region. All error bars correspond to 1-s.d. 
uncertainties. 



the ratio of the sum of the population in the signal and 
idler modes to the reference mode population is below 
10% for small N < 200 and grows to approximately 55% 
for the largest N. In figure [4^i we plot the measured up- 
per bound for the two-mode quadrature variances versus 
local oscillator phase A 2 X± ) ((p) for an exemplary total 
atom number. The sum of the minima at cpo and <^o — § 
provides an upper bound for the sum of the orthogonal 
two-mode quadrature variances EA 2 = A 2 X^ b + A 2 Y+ b , 
while the measurement _of the single- mode quadrature 
variances EA 2 = 2(A 2 X^ + A 2 X^) has been discussed 
above. The generated quantum states fulfill the inequal- 
ity EA 2 , < EA 2 for a wide range of total atom numbers 
(Fig. |4Jd) showing that the produced twin-atom state is 
inseparable - the observed minimal ratio is EA 2 /EA 2 « 
0.2. 

For the squeezed vacuum, which is theoretically expected 
in the non-depleted pump regime, the two-mode quadra- 
ture variances should be reduced below the noise level 
of two coherent modes A 2 X^ h (cp) = 2. For the more 
precisely detectable variance of the quadrature difference 



we find the noise minimum for small total N compara- 
ble to this level, limited by the experimental signal to 
noise ratio. In figure^ we show the dependence of this 
noise minimum on the total atom number N - detection 
noise subtraction is crucial here and without it we find 
A 2 X^ b (ip) ~ 17 as the minimal value. 

In conclusion, we have developed an atomic homo- 
dyne detection method, which allows for the measure- 
ment of quadrature correlations of twin-atom quan- 
tum fields. The continuous variable entangled states 
have been produced in a deterministic manner utiliz- 
ing controlled atomic spin interactions - for small to- 
tal atom numbers the observed populations distributions 
are compatible with the atomic two-mode squeezed vac- 
uum. Einstein-Podolsky-Rosen (EPR) entanglement is 
revealed by A 2 (X t - X i )A 2 (Y t + Y i ) < 1 and we find 
4 ±17 for this value, showing that EPR-entanglement for 
atomic quadratures is within reach. In a future extension 
of our scheme the local oscillator might be split before the 
mixing, enabling precise detection also of the two-mode 
quadrature sum [26] . Spatial separation of the modes can 
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be implemented by employing the Stern- Gerlach tech- 
nique, being an alternative approach to recently reported 
experiments [27] [28] . The regime of large average pair 
population is accessible such that highly non-classical 
quantum states might be generated by controlled quan- 
tum spin evolution [25] . 

We note that in parallel to this work, sub-poissonian 
number fluctuations after spin changing collisions have 
been observed by two other groups [29] [30] and it has 
been shown that the generated quantum correlated states 
are useful for noise interferometry beyond the standard 
quantum limit [30] . 



Supplementary Information 
Experimental sequence 

We routinely produce Bose-Einstein condensates 
(BEC) of Rubidium 87 in the low field seeking (F^mp) = 
(1,-1) hyperfine state with an experimental cycle time 
of approximately 40 s. Before the final evaporation ramp 
in an optical trap we turn up a 1-dimensional optical lat- 
tice slicing the atomic cloud into eight pieces. Since the 
lattice potential does not allow for tunneling between dif- 
ferent sites on the experimental timescale, we start with 
eight independent samples with different mean total atom 
number N per site on which the spin dynamics experi- 
ment is done in parallel. 

First we transfer the atoms to the (1,0) state by a radio 
frequency (rf) tt pulse. A homogeneous magnetic off- 
set field of approximately 9G is applied such that the 
second-order Zeeman shift is sufficient to resolve the dif- 
ferent rf transitions in the F = 1 hyperfine multiplet. 
Afterwards we remove spurious population in the (1, ±1) 
states by a Stern-Gerlach technique. Care is taken to 
adiabatically ramp down the offset field to its final value 
of Bq ~ 1.5 G where we actively stabilize it by a flux- 
gate sensor based feedback loop [2T] , The offset field is 
chosen such that both, state selective two-photon pulses 
between the (1,-1) and the (2, 1) states, as well as power 
broadened symmetric rf coupling between the (2,0) and 
the (2, ±1) states are possible. For the homodyning of 
single- mode quadratures local oscillators in (1, =bl) are 
necessary, which we prepare by transferring a part of the 
population via symmetric rf coupling. The next step is a 
microwave tt pulse from (1,0) to (2,0). 
Spin dynamics is initiated by compensating the second- 
order Zeeman and mean field shift (in total 5 = 162 Hz) 
with microwave dressing of the (2,0) level [23] [24]. The 
dressing field is 98 kHz blue detuned to the (1, 0) (2, 0) 
transition with a resonant microwave Rabi frequency of 
8 kHz. This selective technique has the advantage that 
only spin dynamics between the rnp = 0, ±1 states is 
resonant. In all experiments we allow for 22 ms of spin 
evolution in which the mp = ±1 states become popu- 
lated via spin changing collisions. 

For the measurement of the single-mode quadratures we 
now apply a two-photon tt/2 pulse resonant with the 
(l,=Fl) ^ (2,±1) transitions, where changing the phase 
of these pulses is equivalent to a change of the local os- 
cillator phase. 

In order to realize the three-port beamsplitter scheme, 
i.e. to measure the two-mode quadratures after the spin 
evolution, we apply a single rf field for 60 /is, which cou- 
ples the (2,0) and (2, ±1) states. This corresponds to a 
7r/5 coupling pulse, where the pulse angle is restricted by 
experimental considerations described below. For con- 
trol of the local oscillator phase <p (of the mp = state) 
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FIG. 5. Experimental raw data. Absorption picture 
of atoms in the F = 1 (a) and F — 2 (b) multiplet. The 
counting boxes for regions with atomic signal are shown in 
red, the boxes for technical noise measurement in green. The 
numbers indicate the Zeeman sub state, which is the same in 
each row. Due to the opposite sign of the magnetic moments 
these numbers are inverted between F — 1 and F — 2. The 
size of the image in real space is given. Note that the clouds 
for F — 1 atoms are more extended in vertical direction due 
to the imaging sequence [2T] . 



we switch off the dressing field for a variable time t (be- 
tween ms and 4.5 ms) before the coupling pulse, such 
that (p = 27T 5t. In the case of the measurement of popu- 
lation difference iV_ no coupling is applied after the spin 
changing collisions. 

After the experimental sequences described above the 
population in the different Zeeman states is measured. 
State and site selective imaging is achieved by a combi- 
nation of the Stern- Gerlach technique and a high spatial 
resolution imaging system [21] [31]. The Stern-Gerlach 
pulse, which separates the mp states, is aligned carefully 
with the magnetic offset field such that the following ab- 
sorption imaging is a projective measurement of the pop- 
ulation in the Zeeman states mp denned with the quan- 
tization axis along the direction of the offset field. 



Data analysis 

Imaging is done in the high intensity regime and the 
column density (atom number per pixel) is calculated 
from the raw absorption data following the recipe in [32] . 
Care is taken to assure a linear and well calibrated imag- 
ing system [2TJ [31] . To obtain the atom number per state 
and lattice site we define counting boxes in which the 
column density is integrated. Figure [5] shows a typical 
single run picture with counting boxes indicated. Two 
types of technical noise add to the atomic signal in the 
imaging process. The first is unavoidable photon shot 
noise (V A 2 N psn ^7.0 atoms per box) which we deduce 
from the data in each experimental run [2TJ[3T] . However, 
interference fringes due to the coherence of the imaging 
light make absorption imaging very sensitive to mechan- 
ical vibrations and cause excess noise. Even though we 
discard images which show obviously strong interference 



fringes in the regions with no atomic signal (one picture 
out of 50), the remaining ones still have a finite noise 
contribution due to large scale interference structures. 
These are very hard to avoid, but on an ensemble of pic- 
tures their contribution to the detected variance can be 
well measured. Hence we define counting boxes in re- 
gions without atoms of the same size as those used to 
extract the atom numbers (see Fig. [5| and measure the 
excess variance. It corresponds typically to a noise level 
of VT^N^ * 3.9 atoms per box where we checked for 
the spatial homogeneity of the fringe noise over the rele- 
vant regions of the picture. The technical noise is gaus- 
sian and not correlated with the atomic noise to be mea- 
sured. Therefore we can subtract the independently mea- 
sured noise contributions from the measured atom num- 
ber variances. Note that this subtraction is crucial for 
the reported results since the total background noise level 
corresponds to a variance of A 2 A rt ? ch ^125. Without the 
noise correction we would measure sub-shot-noise fluctu- 
ations only for N > 450 where (JV+) > A 2 A/l ech (see 
figure 2 of the main text). Here it is important to note 
that the uncertainties on the experimental data given in 
this manuscript include the statistical uncertainty due to 
the noise subtraction. As a test for our noise calibra- 
tion we analyzed pictures without atoms in the counting 
boxes and found possible systematic errors smaller than 
the statistical uncertainties. 

The dependence of all presented observables on the to- 
tal atom number N is obtained by binning the data in 
intervals of size 5N = 50. 



Inter-site correlation analysis 

The use of a single site resolved optical lattice in our 
experiments has various advantages. Two of them are 
rather obvious, the boost in statistics due to the parallel 
realization of eight experiments and the increased local 
confinement which is important for the validity of the sin- 
gle mode approximation. Here we point out that inter- 
site correlations between the observables can be used to 
extract phase noise contributions stemming from envi- 
ronmental fluctuations. 

At a finite magnetic field B the energy of the different 
Zeeman states shifts differentially with B which leads to 
phase noise in Ramsey type experiments [21 , in which 
the integrated phase difference is mapped onto an observ- 
able population difference. The homodyning experiments 
reported in this manuscript are of similar type and the 
observed atom number differences are in principle sensi- 
tive to environmentally induced phase noise. Magnetic 
field or microwave phase fluctuations in our experiment 
are mainly low frequency and homogeneous over a spatial 
region much bigger than the size of the entire lattice sys- 
tem. Therefore they result in shot-to-shot fluctuations 
which are correlated between different lattice sites k. 
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FIG. 6. Experimentally observed inter-site correlations, a, Initially coherent state. Strong inter-site correlations 
Cki are detected for an initially coherent state, which are caused by spatially large scale technical fluctuations resulting in 
site correlated phase-noise, b, Spin-changing collisions experiment. After spin dynamics a random phase is observed and the 
inter-site correlations vanish within the statistical uncertainty Cki ~ 0. The striking difference in the correlation signal between 
the two measurements shows that environmental noise can neither explain the random phase of the individual modes nor the 
correlations between signal and idler modes observed after spin-changing collisions. 



Figure [6^i shows the measured inter-site correlations for a 
reference homodyning experiment with the local oscilla- 
tor in the (1,-1) Zeeman state and a coherently prepared 
state in the (2,1) Zeeman state. The inter-site correla- 
tions Cki are calculated for the observed atom number 
differences n k = N±^ k — (1 and 2 label the two Zee- 
man states involved) 



kl 



(n k ni) 
An k Ani 



(1) 



where n k = n k — (n k ) is the atom number difference 
in site k corrected for its mean and An k the standard 
deviation of n k . The experimental sequence is similar to 
the spin-changing collision evolution described above but 
spin dynamics was suppressed by tuning the process off 
resonance. Since the relative phases are in this case ini- 
tially denned by the coherent preparation (corresponding 
to a non- vanishing mean of the single-mode quadrature) , 
finite inter-site correlations \C k i\ > are observed. Their 
strength and sign depend on the difference in the total 
population per site N k —Ni which is due to the N k depen- 
dent mean field shift. In our experiment N k fluctuates 
little from shot to shot. Thus, this population dependent 
energy offset results in a fixed difference in the integrated 
relative phases in different sites which causes correlations 
or anti-correlations depending on its magnitude. Using 
the inter-site correlations to remove these technically in- 
duced phase fluctuations we reduce the observed fluctua- 

An 2 

tions ^ fc -f jv 2 fc ) fr° m approximately 50 to 2.5 times the 
noise limit of two coherent states. This shows that al- 
most all environmental phase noise in our experiment is 
of common mode type and can be detected by the inter- 
site correlations. 

In contrast, for a quantum state with initially vanishing 
coherences such inter-site correlations cannot build up 
such that one expects C k \ ~ 0. In figure [6}d we plot the 
inter-site correlations for the signal mode homodyning 
experiment. Here the (2, 1) Zeeman state was populated 



by spin evolution. No statistically significant correlations 
are observed, excluding that the phase is randomized by 
environmental noise. We also tested for inter-site corre- 
lations in the signal-idler mode entanglement measure- 
ments (Fig. 3a, b and 4 of the main text) and found a 
vanishing signal C k \ « excluding a significant influence 
of environmental noise. 

We emphasize that this characterization of classical de- 
phasing is crucial for the results presented in figure 4b of 
the main text, since we use different local oscillators to 
deduce the single- and two-mode quadratures. Stronger 
classical de-phasing in the single-mode quadrature mea- 
surement compared to the two-mode measurement would 
result in a false bound for two-mode inseparability. 

Two- mode quadratures 

The two-mode quadratures of signal (a^) and idler (a^) 
modes are measured by simultaneous symmetric radio- 
frequency coupling to the pump mode ao which serves as 
the local oscillator. This process is described by a uni- 
tary rotation U cp \ = Q~' lH ^ t / h i n the three mode system 
generated by the Hamiltonian 

# cp i = ^^( a °4 + a o a t + a l a i + a o4) ( 2 ) 

where 2ttH is Planck's constant. 

Experimentally the coupling is switched on and off 
non-adiabatically with a total pulse duration r cp i . The 
rotation (mixing) angle of the pulse is given by 7 = Q r cp i , 
which is defined such that a 7 = tt pulse transfers all pop- 
ulation from the mp = state to the mp = ±1 states 
(assuming they are initially empty). 
We choose a mixing angle of 7 « 7r/5, since we have to 
make a tradeoff between signal-to- noise ratio (optimal for 
7 = 7r) and an unwanted spurious coupling between the 
mf = ±1 and = ±2 states, favoring a smaller mixing 
angle 7. 
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The unitary transformation generated by the coupling 
Hamiltonian ^ is described by the matrix 

ov I (3) 



U 



( 2 



c 7 + l 



72 



with s 7 = sin(^), c 7 = cos(^). The mode oper- 
ators (a^,a ,a^) T in the Heisenberg-picture after the 
coupling pulse are obtained by applying the matrix U 
to (a^,ao,a^) T . The phase shift <p of the local oscil- 
lator before the pulse is described by the substitution 
ao —> ao((p) = a^e lip . For the sake of clarity this abbrevi- 
ation is used in the following without the explicit phase 
dependence. 



Two-mode quadrature difference 



For the experimentally measured population difference 

NL{<p) 



— " / ^ ,/ — a'^a^ in the m F = ±1 states we find 



N - = a i a i 



c 7 7V_ 



-|- (2^(2o — aj^ao 



a t4) 



(4) 



°7 ( ~t 

The second term is, up to a | phase shift of the ao mode, 
proportional to the generalized two-mode quadratures 



X_((p) = (a+aj - 



a\a 



• h.c.)/ V (oja ) 



(5) 



In the limit of a large local oscillator (ao ~ \/ (No)e l(p ) 
the generalized two-mode quadratures correspond to the 
canonical ones [26]. With the normalization factor 
(iVo)s^/2 and by using (iV_) = and t 7 = tan(^) we 
obtain the measured upper bound 

A 2 X» h (ip) = A 2 N f _&)/((N )s 2 7 /2) 

= A 2 X_((^) 

2 9 
A N- 



(6) 



:(iV_X_(^)+X_(^)iV_) 



The symbol A 2 indicates the variance of the subse- 
quent variable. The local oscillator phase independent 
variance offset (second term) which is proportional 
to the variance of the population difference A 2 7V_ is 
small (approximately 0.1) for small total atom numbers 
(150 < N < 200) and increases to approximately 3.2 for 
500 < N < 550. The terms in the last line vanish due 
to their strong sensitivity to magnetic field fluctuations 
(700Hz/mG). At this sensitivity the coherence time in 
our experiment is approximately 1.5 ms which results in 
(N-X-) w after 22 ms of spin evolution. 



Two-mode quadrature sum 



The second important two-mode quadrature 



(7) 



i.e., the sum of the two single mode quadratures, can be 
measured less precise in the current scheme, however, an 
upper bound for its variance can be obtained. For the 
sum of the population in the mp = ±1 states = 



a} a\ + a'laL we obtain 



,'t„/ 



N+ + 

*- l < t 



— a^ao 



ajao 



(8) 



The third term zc j| 7 (a^a^ — ajao — ajao + cl^Oq) is pro- 
portional to the two- mode quadrature X+(p-\- 1), which 
is the only local oscillator phase dependent term. Covari- 
ance terms between this one and the three other terms 
in equation ([§]) vanish due to their strong sensitivity to 
magnetic field fluctuations (see above). Hence, one ob- 
tains the upper bound for the two-mode quadrature sum 



A 2 Xf {ip) = A 2 N' + (p)/((N )^ s 2/2) 
= A 2 X + (p) + const. 



(9) 



where the constant offset due to the variance of the phase 
independent terms in equation ^ limits the obtainable 
precision when deducing A 2 X+(cp) < A 2 N^_(cp) from our 
measurements. Note that the normalisations for X + and 
X_ are different. A future extension of this scheme to 
measure also the sum of the single mode quadratures with 
high precision and therefore both quadrature Einstein- 
Podolsky-Rosen variables [7; would require the splitting 
of the local oscillator prior to the coupling, such that the 
phase of the two coupling fields can be tuned indepen- 
dently [33] . 
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